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We study a generalisation of the Hatano-Nelson Hamiltonian in which a periodic modulation of 
the site energies is present in addition to the usual random distribution. The system can then 
become localized by disorder or develop a band gap, and the eigenspectrum shows a wide variety 
of topologies. We determine the phase diagram, and perform a finite size scaling analysis of the 
localization transition. 

PACS numbers: 71.10.Fd, 71.30.+h 



I. INTRODUCTION 

The Anderson transitioni'^, in which eigenstates of 
non-interacting electrons are locahzed by disorder, has 
long been a paradigm of a metal-insulator transition in 
solid state physics. The metallic phase occurs only in 
three dimensions, since any amount of randomness lo- 
calizes all particles in one or two dimensions. In three 
dimensions, a mobility edge separates localized states 
whose energies lie at the edges of the spectrum, from 
extended states whose energies lie in the middle. The 
mobility edges appear for a critical value of the disorder 
and move with the disorder until no extended states re- 
main. If the Fermi energy is located below or above the 
mobility edge, inside the localized states, the system is 
then an insulator. Recently, research on this subject has 
experienced renewed interest as localization of cold atoms 
in the presence of a disordered optical potential has been 
observed^ii. Whereas experiments in solid state physics 
are intrinsically limited as interactions between electrons 
cannot be neglected, cold atoms allow the exploration of 
the non-interacting regimes with possibilities of control 
that could not be achieved before. 

Another, simpler origin for an insulating behavior in 
a solid is the presence of gaps due to a periodic poten- 
tial applied to the particles. These two possible insulat- 
ing phases are in competition with each other as they 
are based on opposite behaviors. The Anderson insula- 
tor is due to the localized nature of the wave function, 
whereas the gapped insulator is due to its extended na- 
ture. A third source of insulating behavior is interactions 
between particles. In many ways this is the most difficult 
situation, since it involves treating the many-body prob- 
lem accurately. In the systems with both interactions 
and randomness, competition between Mott and band 
insulating phases have been shown to result in metallic 
phases, as the two sources of insulating behavior coun- 
teract each otheri^. The possibility of metallic phases 
arising due to the addition of interactions to two dimen- 



sional disordered systems has also been an area of great 
intereslj^Ti^. 

The goal of this article is to study the interplay 
between the two sources of single particle insulating 
behavior- disorder and gaps due to a periodic potential- 
in the context of non-Hermitian matrices. For this we 
will use a modified version of the Hatano-Nelson model 
(HNM), that includes at the same time a periodic po- 
tential that generates gaps in the eigenspectrum and a 
disordered potential that will drive a localization tran- 
sition. The HNMiSiii is a single particle, random, non- 
Hermitian Hamiltonian that was introduced to study the 
motion of magnetic flux lines in disordered type II super- 
conductors, the path integral representation of the HNM 
being analogous to a classical model used to describe the 
flux lines. The HNM shows a localization transition in 
one dimensioniSrJ^. It is then a very convenient tool to 
study the localization transition numerically, as it allows 
the study of systems with large linear sizes. 

The paper is organised as follows. In the Sec. II we 
introduce the HNM and recall some basic properties. In 
Sec. HI, we present the different regimes which can be 
observed. Sec. IV contains a detailed analysis of some of 
the transitions appearing in the generalized HNM. The 
conclusion describes some of the connections the addi- 
tional term we introduce into the Hatano-Nelson model 
has with the eigenspectrum of several physical systems 
for whose properties random matrices have been used. 



II. MODIFIED HATANO-NELSON MODEL 

A. Hatano-Nelson model 

The Hatano-Nelson model is a discrete model similar 
to the Anderson model used to study localization. In one 
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Here cj,(c^) are the usual fermionic creation (destruction) 
operators and the number of particles on 

site X. The first term of the Hamiltonian describes the 
hopping of particles between sites, and the parameter h 
controls the asymmetry between the hopping amplitudes 
in the left and right directions. We will choose the ran- 
dom sites energies fix according to a uniform distribution 
on the interval [—A/2, +A/2]. The hopping parameter t 
is set to one to fix the energy scale and L is the number 
of lattice sites. 

In the absence of disorder (A = 0) the eigenenergies 
are given by, 
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where k = 1,2 •••L. All the eigenvalues (except k = 
i/2, L) are complex and lie on an ellipse in the complex 
plane, centered at the origin. When disorder is turned 
on, A 7^ 0, the eigenvalues of the d = 1 HNM remain con- 
fined to one dimensional curves in the complex planei^, 
however, some of the eigenvalues become real. The corre- 
sponding eigenstates become localized (Fig. [T]) . We will 
continue to use the notation for the eigenvalues for 
nonzero A even though momentum is no longer a good 
quantum number due to the breaking of translation in- 
variance by the disorder, k will be understood to be a 
generic mode label. 
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One way to understand the appearance of real eigen- 
values is to consider a system with open boundary condi- 
tions (OBC)iS'ii. In this case, all the eigenvalues of 'Hhn 
are real, despite the fact that the Hamiltonian is not Her- 
mitian. In fact, the eigenspectrum is independent of /i, 
since the e^'* factors can be "gauged away" by a suitable 
redefinition (f){x) = e~''^0(a;), of the eigenvector compo- 
nents. The HNM spectrum becomes identical to that of 
the Anderson model. On a system with periodic bound- 
ary conditions (PBC), however, this eigenvector trans- 
formation is no longer possible, and some eigenvalues are 
complex (see Fig. [T]). 

The dramatic difference in the spectrum between OBC 
and PBC of course reflects a rather general property of 
the localization transition: A localized state is insensitive 
to the boundary conditions and its eigenvalue should re- 
main real, even in the presence of PBC. On the contrary, 
an extended state, for example a plane wave state, is 
sensitive to the boundary conditions and its eigenvalue 
should be complex in the PBC case (Fig. [T]). 

Le Doussali^ has further clarified the relation between 
the eigenvector tranformation and the localization length 
by pointing out that the eigenvalues will remain real if 
the localization decay e~'/^ is more rapid than the growth 
e''' induced by the transformation. As a consequence, the 
localization length ^ does not diverge at the transition, 
but takes on the value ^ = l//i at the bifurcation point 
in Fig. [1] where the line of real eigenvalues splits off into 
the complex plane. 

For the particular case of a Cauchy distribution of 
randomness, P{A) = 7/7r(7^ + A^) the location of the 
bifurcation point is knowni^ to be given by Re[efc] = 
•\/ {2t sinh h)^ — j'^ / tanh h, and Refe^-] = at the crit- 
ical value 7c ~ 2tsinh/i. At this point the randomness 
is sufficiently large that all eigenvalues are real, and all 
modes are localized. 

More direct evidence of the localized nature of the 
eigenstates with real eigenvalues can be obtained by di- 
rectly looking at how the corresponding density varies 
with distancei^. For example, we can study the partici- 
pation ratio Pfc for a given eigenvalue e^. As the Hatano- 
Nelson model is non-Hermitian, the left eigenvector (fc| ~ 
J2x(^\'l^k(^) and right eigenvector \k) = J2x '/'fc (^)N) 
not Hermitian conjugates. It was showni^ that the prob- 
ability density in these systems is given by (a;)0^(a;)| 
which is a better behaved quantity than |0^(a^)P or 
|(/)^'(x)p although they give quite similar results for lo- 
calized states (Fig. [3]). 

Then the participation ratio can be defined as. 
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FIG. 1: (Color online) Typical spectra of the Nelson-Hatano 
model. In the presence of disorder, localized states (real eigen- 
values) appear for A larger than a critical value. Half of the 
spectra are shown and the spectra are approximately sym- 
metric around the imaginary axis. The two larger up triangles 
correspond to the eigenstates used for Fig. [2] 



Pfc corresponds to the portion of the lattice sites that is 
occupied for a given eigenvalue. It goes to zero for a lo- 
calized state in a large system and goes to a number of 
order unity for an extended state when the density dis- 
tribution is spread equally over all the sites. Fig. [3] shows 
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FIG. 2: (Color online) Three different possible definitions of 
the density in an extended (left) and localized (right) states 
corresponding to different eigenvalues of the same sys- 
tem. In general \(j>k'^i^)\^ show quite strong fluctuations, 
thus making \(f>^{x)cj>^{x)\ a better quantity to identify the 
different states. 



FIG. 3: (Color online) Different quantities which distinguish 
the extended from the localised states in the Hatano Nelson 
model: the imaginary part of the eigenvalues Im(efe), partici- 
pation ratio Pfc, and real Re(Jfc) and imaginary Im(Ji;) parts 
of the current. All these quantities give consistent information 
concerning which modes are localized or not. 



that the vanishing of the participation ratio aligns nicely 
with the energy at which the eigenvalues are real. Similar 
participation ratios can be defined with other definitions 



of the density P^''^ = ^/L. But, as for the 

densities (Fig. ^ themselves, Pk is better behaved than 



The current Jk defined as Ju = —idek/dh^ is com- 
pared with the participation ratio in Fig. [31 It gives sim- 
ilar results for the extended or localized nature of the 
states, but is more difficult to calculate as it requires a 
numerical derivative (and then two diagonalization in- 
stead of one for the participation ratios) and is difficult 
to define close to the mobility edges. 



B. Effect of a periodic potential 

Our interest here is in adding a simple potential of 
period 2 to the original HNM, 
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FIG. 4: (Color online) With the periodic potential, a gap is 
opened for V > V° — tsinh(/i) although all the states remain 
extended. ~ 0.636 in the case shown, t — 1.0, h = 0.6. 



When h = A = 0, this term generates a spectrum with 
two bands separated by a gap of width 2V. fn the non 
disordered case where h and A = we obtain the 
following eigenvalues 



given by 



Wc 



2JV^-f^ sinh2(/i) 
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with k = 1,2 - ■ ■ L/2. In this ordered case, the spectrum 
has a gap if V > V° = isinh(/i) and its width Wg is 



Even when A = 0, there is a transition in this model 
corresponding to the opening of the gap at half-filling 
(see Fig. HI . The gap opens with exponent /3 = 1/2 since 
near the transition Wg ~ yy2V°{V - V°). 
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FIG. 5: (Color online) Different topologies of the eigenspec- 
trum for one realization of a L = 1000 system and for h = 0.6. 
The simple E topology (only extended states without a gap) 
and the rarely observed LELEL are not shown. See text for 
discussion. 



III. PHASE DIAGRAM 

In the presence of both disorder and the periodic po- 
tential, the spectrum of the model can show several dif- 
ferent topologies. We calculate the eigenspectrum for one 
disorder realisation of a large size system (L = 1000) and 
analyse the different topologies appearing. We checked 
that these topologies were still present for other sizes 
(L = 500, i — 2000) and other disorder realisations. 

Scanning horizontally along the real axis in the com- 
plex eigenvalue plane, we observe three types of regions 
in the spectrum. First we observe cases of real eigen- 
values corresponding to localized (L) states. Second we 
observe regions with complex eigenvalues corresponding 
to extended (E) states. Finally we observe some gaps (G) 
in the spectrum. A given topology is characterized by the 
succession of these different regions as one increases the 
real part of the eigenvalue, Kc{ek). The different topolo- 
gies we observed are shown in Fig. [S] and their position 
in the {V, A) plane for a fixed value h = 0.6 is shown in 
Fig. [51 As an example, the phase shown in Fig.[5l^a) is de- 
noted as the LEL phase since we observe a central region 
of extended states (E) between two regions of localized 
states (L) with no gaps. 

The location of these different phases for h — 0.6 is 
shown in Fig. [5] for a given realization. As expected, a 
gap appears for V large enough. For moderate disorder 
(A < 2), the opening of the gap happens approximately 
at the same value V° = sinhO.6 = 0.636 as in the or- 
dered, A = case. Upon increasing A, the original low 
disorder E and EGE phases are modified as some of the 
extended states are transformed into localized ones at 
the edges of the extended regions, leading to the LEL 



FIG. 6: (Color online) The different phases observed in the 
Hatano-Nelson model for a given realization on a L = 1000 
system size and for t — l,h = 0.6. At the boundary between 
different phases, the behavior is sometimes irregular for a fi- 
nite size system and does not fall into one of the identified 
categories. The labels (a)-(f) correspond to the cases shown 
in Fig. El 



and LELGLEL (or LEGEL) phases respectively. In the 
LEGEL phase, localization first appears on the external 
edge of the extended regions and not on the border of the 
gap. This is consistent with the fact that localized states 
appear first for the extremal eigenvalues in the Anderson 
model. 

In a very small region of the phase space, we also ob- 
serve a LELEL phase where the region between the two 
extended blobs is filled by localized states but no gap 
appears. As A is increased further, all the extended 
states disappear leading to the large disorder L and LGL 
phases. When V and A are both large, the transition 
from the L phase to the LGL phase occurs at Vc = A/2 
which is the result expected in the one site limit, when 
the hopping terms are neglected. In the one site limit 
the gap opens like Wg — 2(T^ — T4) which is a different 
behavior from the opening of the gap due to the periodic 
potential. 

For a given realisation, close to the boundary between 
two phases, the spectrum becomes disordered and it is 
difficult to classify it as belonging to one of the simple 
topologies identified before. These regions are indicated 
by "x" in the chosen example (Fig.|6l). To calculate more 
precisely the boundaries between the different phases, it 
is then necessary to use different quantitative indicators 
and to average those over several disorder realisations. 
First we measure the width of the gap in the system and 
determine whether it is present or not. Second we deter- 
mine if the spectrum is made of complex values (extended 
states) only, of a mixture of complex and real values (ex- 
tended and localized states), or of real values only. We 
consider a L = 1000 system and average over 10 real- 
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FIG. 7: (Color online) The phase diagram obtained by aver- 
aging several disorder realisations to calculate the boundary 
between different regions. The smooth lines are fits to the 
data. The continuous green line marks the limit between 
gapped and ungapped regions. The dotted blue, the limit 
between a spectrum without any real values and a spectrum 
with real and complex values. The dashed red line is the limit 
between a spectrum with real and complex eigenvalues and a 
spectrum with only real values. 




FIG. 8: (Color online) For V = 0.7, upon increasing the disor- 
der, the gap Wg is closed before the states becomes localised, 
which is shown by the participation ratio P becoming zero. 
Here, the quantities are averaged over only 5% of the total 
number of eigenstates located in the middle of the spectrum, 
to get rid of possible bias induced by states located at the 
edges of the spectrum. The measurements show the cross- 
ings of the two boundaries of the phase diagram at largest A 
(gap closing and all extended states vanishing) for V = 0.7 
in Fig. [T] but not the crossing of the boundary at smallest A 
when localized states first appear. 



isations. Drawing these three lines in the phase space, 
we obtain the phase diagram shown in Fig. [71 With 
these quantities, we cannot distinguish the LELGLEL 
from the LEGEL phases, as weh as the LELEL from the 
LEL phase. 

We observe an interesting reentrant behavior. If the 
system is in the EGE phase but the gap is not too large, 
when we increase A, the gap will be closed by the dis- 
order but the states in the middle of the spectrum will 
remain extended (Fig. |5]). Then we observe a transition 
where increasing disorder drives the gapped system into 
an ungapped phase. This can be understood as follows: 
as the disorder is increased, the E regions have a tendency 
to flatten out as they evolves towards L regions with only 
real values. By doing so, the E regions spread along the 
real axis. If the gap is not to large, two gapped E regions 
can then merge in the process, which is happening here. 



transition between those two different behaviors happens 
precisely at the point where, in the phase diagram Fig. [71 
we are passing from the low disorder into the large disor- 
der regime. For A < 2, the transition almost takes place 
at the non- interacting critical value V°, the transition 
occurring at nearly constant V and with the square root 
behavior. Above this value, the transition takes place ap- 
proximately at Vc = A/2 with the linear behavior, which 
is characteristic of the disordered regime. 



B. Finite size analysis of the localization transition 



IV. TRANSITIONS 
A. Opening of the gap 

As expected, we observe two different behaviors for 
the opening of the gap for a fixed value of A as 1^ is 
increased (see Fig. U 6 ^ V — Vc, we have a behav- 
ior Wq (X VS as in the non interacting regime, whereas 
Wg (x S in the strongly disordered regime Fig. [SI The 



It is possible to analyze more precisely the transition to 
the localized regime by performing a finite size analysis of 
the participation ratio. We define the mean participation 
ratio P = Sfe-Pfc/^ which becomes zero when we have 
only localized states. We found that the participation 
ratio is varying linearly with the size of the system L. In 
Fig. [TUl we show the result of a finite size analysis for a 
given value of V and for different values of A. We see that 
P extrapolates to zero for values larger than Ac ~ 2.5. 
This is in agreement with the phase diagram obtained 
with a L = 1000 size (Fig. [7]). 
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FIG. 9: (Color online) Evolution of the gap as a function of 
V. For V < 2.0, the gap opens like Wq oc y/V — K as in the 
non disordered case. For V > 2.0, the opening of the gap is 
linear as expected for a transition in the disorder regime. 
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FIG. 10: (Color online) Finite size scaling analysis of the 
participation ratio. For this value of V, the transition towards 
localization happens at Ac — 2.5. 



V. CONCLUSION 

We have studied the Hatano-Nelson model in the pres- 
ence of a periodic potential. We have determined the 
phase diagram of the system and shown that it presents 
a variety of different localized or gapped phases, distin- 
guished by the topology of the eigenspectrum in the com- 
plex plane. Interesting insulator-metal transitions can 
occur, driven by increasing disorder. Finally, we have 
determined how the gap and participation ratio evolve 
close to these transitions. 

There are some interesting analogies between eigen- 



A = with hoppings , ty and non-Hermiticity param- 
y. In fact, the eigenvalue distribution in the 



spectra features induced by a periodic potential, dis- 
cussed here, and those which appear in crossing over 
from one to two dimensionsiS. We comment first that 
the d = 2 model can be trivially solved in the clean limit 

eters h^, h 

complex plane is obtained by taking ellipses determined 
by ty and hy, and placing them with centers at each of 
the eigenvalues of Eq. [2] In the limit of large lattice 
sizes these eigenvalues become space filling in the com- 
plex plane. 

However, just as discussed here, different topologies 
are possible. If ty and hy are small compared to and 
hx (or vice-versa) one is overlaying small ellipses on a 
background of a larger ellipse, and the eigenspectrum re- 
tains a hole in the origin. When the hoppings in the two 
dimensions are comparable, in contrast, the central hole 
can get filled in. It is not surprising that a particularly 
close analogy is obtained by considering a system consist- 
ing of two chains since such a situation describes, as with 
the periodic potential considered here, a lattice with two 
different types of sites (two bands). The resulting eigen- 
value spectra reflect considerable similarities with those 
found in this paper. See, for example. Fig. 2 of [l6| . 

Random non-Hcrmitian matrices also arise in the 
study of quenched Quantum Chromo-Dynamics at finite 
quark densityiiiii^ While the structure of the matrices 
does not take the precise form of the HNM, it is inter- 
esting that in certain cases ( "Dyson index" /3 = 4) a gap 
appears in the spectrum in the complex plane, analogous 
to that found hereJ^ That is, the eigenvalue distribu- 
tion can be driven to zero in a region in the vicinity of 
the imaginary axis, not unlike the effect of increasing V 
in Fig. m (In these studies the spatial dimensionality is 
greater than one so that the eigenvalues are space fill- 
ing as opposed to lying along one-dimensional curves.) 
Indeed, a key result^i is that the appearance of a gap 
around the imaginary axis requires a finite value of the 
appropriate tuning parameter, the chemical potential /i, 
much as our phase diagram (Fig. requires a finite stag- 
gered potential V > Vc — 0.7 for the appearance of a gap. 

The original motivation for the HNM was the problem 
of flux line motion in type II superconductors. Specifi- 
cally, deliberate heavy ion irradiatior>i^ produces random 
tracks along which the flux lines like to sit, and a trans- 
verse magnetic field then provides the 'non-Hermiticity' 
which attempts to delocalize them. It is also possible to 
provide regularly patterned pinning centers, which would 
be a realization of the staggered ±V term in the exten- 
sion of the HNM presented here<^"— Our phase diagram. 
Fig. ini indicates that the amount of random disorder A 
required to produce localized phases in general increases 
with the amplitude V of the patterned potential, so that 
critical currents would decrease as patterned pinning cen- 
ters are established if the modulation is at the atomic 
level. It would be interesting to study the effect of reg- 
ular variation in the site energy at longer length scales, 
such as would more typically be accessed in experiments 
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with patterned pinning. 
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